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Abstract 



This paper is a continuation of our paper [8J in which we proved the Central Limit Theorem 
for the matrix elements of differential functions of the real symmetric random Gaussian matrices 
. (GOE). Here we consider the real symmetric random Wigner matrices having independent 

(modulo symmetry conditions) but not necessarily Gaussian entries. We show that in this case 
the matrix elements of sufficiently smooth functions of these random matrices have in general 
another limiting law which coincides essentially with the probability law of matrix entries. 

Oh 

! 1 Introduction 

> 

& . 

We are interested in asymptotic properties of matrix elements ipj k (M), j,k = l,..,n, n — > oo, 
where ip is a smooth enough test-function, and M is the Wigner matrix. We define the Wigner real 
■ symmetric matrix as follows: 

|>; M n = n-^W n , W n = {W$eR,W$=Wlf}% k=1 , (f.l) 

m : 

CN ■ where {W^ }i<j<k<n are independent random variables satisfying 

cn 

E{W$} = 0, E{(W$) 2 } = w 2 (l + 6 jk ). (1.2) 



The case of the Gaussian random variables obeying (|1.2p corresponds to the Gaussian Orthogonal 
Ensemble (GOE) (see e.g. [5]): 



M n = n-^Wn, W n = { W jk = W kj G R, W jk € Af(0, w 2 (l + ^))}™ fc=1 . (1.3) 

(n) 

We will assume in what follows additional conditions on distributions of Wj k , mostly in the form 

(n) 

of existence of certain moments of W), , whose order will depend on the problem under study. 

In our paper [8j we have considered matrix elements of functions of the GOE matrices and have 
proved the following facts. 

Theorem 1.1 Let M n be the GOE matrix \1.3\) , and ipi2 : R — >• C be bounded test functions with 
bounded derivative. Denote 

{^{M n ))° n = (<Pi(K))jj ~ E{(<Pi(K))jj} (1-4) 

and 

Cov{(^i(M„)) ji , M M n))jj} = E{(^i(M n )) J - i ( ¥?2 (M n ))° i }. 
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Then we have for any j n = 1, n 

limnCovj^!^)) ,(<p 2 (Mn)) } 

n— >00 V 'Jnjn \ ' Jnjn J 

/ A^iA^ 2 p sc (Ai) / o sc (A2)dAidA 2 , (1.5) 

-2w J-2w 
where 

Acp = (f(Xi) — ¥>(A 2 ), (1.6) 

/> sc is t/ie density of the semicircle law 

p sc (\) = (2Trw 2 )- 1 (4w 2 -\ 2 ) 1 / 2 , (1.7) 

and x+ = max{0,x}. 

Theorem 1.2 Let M n 6e £/ie GOE matrix \1.3\) . and ip : R — >■ R 6e bounded function with bounded 
derivative. Then for any j n = l,...,n the random variable y/nip°(M n )j n j n converges in distribution 
to the Gaussian random variable with zero mean and the variance 

f -2w f -2w 

tGOE\ 



Vd w\ = / / (A^)^ sc (Ai)p sc (A 2 )dAidA 2 . (1.8) 

In the present paper we prove a counterpart of Theorems 11.11 and 11.21 for the Wigner matrices. 
In particular, we show in Theorem 13.21 below that in this case the r.h.s. of (|1.5p has an additional 
term proportional to the fourth cumulant K/± of non-diagonal entries (see (|2.18p for the definition). 
This result is in accordance with that for centered linear eigenvalue statistics 

n 

KW] = TiV(M„) = X)^(M n ) (1.9) 
i=i 

of the Wigner matrices (see Theorem 3.6 of [?])• On the other hand the individual matrix elements 
y/n<p'jj(M) do not satisfy in general the Central Limit Theorem (CLT). In Theorem !3.4l of this paper 
we find limiting probability law for yJnyfjAM) which is not Gaussian in general but rather that of 
the sum of the probability law of entries of M n modulo a certain rescaling. To obtain the CLT, one 
has to impose an integral condition on the test function, i.e., the set of test functions for which we 
have the CLT has the codimension one. 

Our results of [8j and of this paper can be viewed as analogs of the E. Borel theorem on the 
limiting probability law of entries of orthogonal matrices of size n as n — > oo (see e.g. [2]). 

Convention: We will use letter C for an absolute constant that does not depend on j, k, and n, 
and may be distinct on different occasions. 



2 Technical Means 

To make the paper self-consistent, we present here several technical facts that will be often used 
below. For the proof of these facts see e.g. [U [7] 

We start from the generalized Fourier transform, in fact the tt/2 rotated Laplace transform (see 
e.g. [15] . Sections 1.8-9 for its definition). 

Proposition 2.1 Let f : R + — > C be a locally Lipshitzian and such that for some 5 > 

supe-**|/(<)l <oo, (2.1) 
t>o 
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and let f : {z £ C : 9z < —5} — > C be its generalized Fourier transform 

____ poo 

f(z) = r 1 / e~ lzt f(t)dt. (2.2) 
J o 

TTie inversion formula is given by 

f(t) = ±je izt f{z)dz,t>0, (2.3) 

where L = (— oo — ie, oo — e > <5, and the principal value of the integral at infinity is used. 

Denote for the moment the correspondence between functions and their generalized Fourier trans- 
forms as f /. Then we have: 

(i) tif(T)dT^(iz)- l f( Z )- 

(ii) /„' h(t - r)f 2 (r)dr := (A * f 2 )(t) o t/i(z)/ 2 (z); 
(mj «/-P, Q, and -R are differentiable, and R(0) = 0, then the equation 

P(t) + I dh I' Q(h - t 2 )P(t 2 )dt 2 = R(t), t > 0, (2.4) 
Jo Jo 

has a unique differentiable solution 

P(t) = - [ t T(t-t 1 )R'(t 1 )dt 1 , (2.5) 
Jo 

where 

T^iz + Q)- 1 (2.6) 

provided by 

z + Q{z) / 0, ^sz < 0. (2.7) 
The next proposition presents simple facts of linear algebra 

Proposition 2.2 Let M and M' be n x n matrices and i £ R. Then we have the following: 
(i) the Duhamel formula 

e (M+M')t = e Mt + /"* e M(t-s) M , e (M+M>)s ds _ (2 _ g) 
JO 

(ii) if for a real symmetric M 

U(t) = e itM , t el, (2.9) 
then U (t) is a symmetric unitary matrix and 

n 

U(t 1 )U(t 2 ) = U(h + t 2 ), ||tf(t)|| = l, J2\U jk (t)\ 2 = 1, (2.10) 

i=i 

so that 

n n 

\u jk (t)\ < i, ^ nV2 ' IMW'^)! < i; (2.ii) 

fc=i fc=i 
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(Hi) 

D jk U ab {t) = i/3 jk (U aj * U bk + U bj * U ak ) (t), (2.12) 

where 

D jk = 8/dM jk , (2.13) 

Pjh = (1 + ^'fc)" 1 = l-Sjk/2, (2.14) 
the symbol "*" is defined in Proposition ^. 1\ (ii), and 

\D l jk U ab (t)\ < ci\t\ l , Cl = 2 l /l\. (2.15) 

Now a generalization of a property of the Gaussian random variable £ of zero mean and of 
variance w 2 according to which if <I> : R — > C is a differentiable function with polynomially bounded 
derivative, then 

E{£$(0} = E{£ 2 }E{$'}. (2.16) 

Formula (|2.16p is a particular case of more general formula. To write it we recall some definitions. 
If a random variable £ has a finite pth absolute moment, p > 1, then we have the expansions 

f(t) :=E{e^ = f| !t y + o(f), 
j=0 J ' 

and 

/(£):= log E{e^} = £^) J +°(n t -+ 0, (2.17) 



where " log " denotes the principal branch of logarithm, the coefficients in the expansion of / are 
the moments {Mj} of £, and the coefficients in the expansion of I are the cumulants {kj} of £. For 
small j one easily expresses Kj via //i, /j>2, ■ ■ ■ , [J-j- In particular, 



«l = Ml> k 2 = n 2 ~ Mi = Var{£}, k 3 = /i 3 - 3^>Mi + 2 Mi\ ( 2 -18) 

In general 



K 4 = M4 - 3/U2 - 4/i 3 ^i + 12/^Mi - 6yu|, ... 



%=^c A ^ A , (2.19) 



A 

where the sum is over all additive partitions A of the set {1, . . . c\ are known coefficients and 
MA = Wiex^h see e.g. [H]. We have 

Proposition 2.3 Let £ be a random variable such that E{|£| p+2 } < 00 for a certain non-negative 
integer p. Then for any function <E> : R — > C of the class C p+1 with bounded partial derivatives $0, 
/ = 1, + 1, u;e Ziaue 

Etf*(0} = E^f E{^ (/) (e)}+^ (2.20) 

w/iere 

\s p \ < C p B{\e +2 }^P\^ P+1) (t)\, C p < l + f + /f7 +2 - ( 2 - 21 ) 

If the characteristic function E{e J *l^} is whole, and <I> £ C°°, t/ien 

00 

E{^(o} = E^r E{$(0(e)} (2 - 22) 

«=0 
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provided that for some a > 

|E{$«(£)}|<a z , (2.23) 
and for some R = ca, c > 1, K\ satisfy the condition: 

£k^<oc. {2M) 
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1=0 

3 Covariance of Matrix Elements 

We show first that if M is the Wigner matrix with uniformly bounded sixth moments of its entries, 
and the test-function 93 is essentially of class C 3 , then the variance of y/nip^{M) is of the order 
O(l) as n — > 00. We have 

Lemma 3.1 Let M = n~ l / 2 W be the real symmetric Wigner matrix M-l\) - Al.fy) . Assume that: 

(i) the third moments of its entries do not depend on j, k, and n: 

n 3 = E{(W^) 3 }; (3.1) 

(ii) the sixth moments are uniformly bounded: 

w 6 :=sup max E{(W% ] f} < 00. (3.2) 

nSN l <h k < n 

Then for any test-function ip : M — > C, whose Fourier transform 

1 I" -itX, 



m = ^J e- lt V(A) d\ (3.3) 
satisfies the condition 

J(l + \t\f\fi(t)\dt<oo, (3.4) 

we have the bound 

Var{V^-(M)} : = E{|V^-(M)| 2 } 

<c(J(l + \t\f\(p(t)\dty. (3.5) 

Proof. It follows from the Fourier inversion formula 

i\t : 



<p(\) = j e lM (p(t)dt (3.6) 
and the spectral theorem for symmetric matrices that 

^•(M) = / £(t)U° 3 (t)dt, (3.7) 



where U is defined in ([2.9p . This and the Schwarz inequality yield 

Var{^(M)} = j j ^{UjjWUJjWffltlWWdtidh 
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Now (13.51) follows from the estimate 



'33 

proved in Lemma [3.81 below (see Appendix, (|3.106p ). and condition f|3.4j) . 



(3i 



Var{Ujj{t)} < C(l + \t\f/n (3.9) 
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Theorem 3.2 Let M = n l / 2 W be the real symmetric Wigner matrix - (l.ty) . Assume that 

the third and fourth moments do not depend on j, k, and n: 

M3 = E{(^) 3 }, ^4 = E{(^) 4 }, (3.10) 

and the sixth absolute moments are uniformly bounded (see \3.2\) ). Let <pi >2 : K — > C be the test- 
functions, whose Fourier transforms pi >2 fS satisfy flff.^D - Then we have for any j = j n £ [l,n\: 

/2w p2w 
/ A^iA^ 2 p sc (Ai) / o sc (A 2 )dAiciA2 
-2w J-2w 

2 r-2w 



-in w(\)(w 2 - \ 2 ) Psc (\)d\, (3.H) 
w i=iJ- 2w 



where A(p is defined in M.6\) , and 

K4 = ha — 3w 4 (3-12) 

is the fourth cumulant of the off-diagonal entries (see A2.18\) ). 
In particular, 



V^[(p] : = lim Var{v^(M)} 

= V d GOE [p] + ^ r V(»)(w 2 ~ f)PsMdn (3.13) 



2w 



withV^ Ub [p\ of fPj) . 

Remark 3.3 (%). If p is odd, then V d w \<p) = V<f OE [cp}. 

(ii). Note that we choose here the Wigner matrix so that its first two moments matches the 
first two moments of the GOE matrix (see ([1.2p ). This fact allows to use known properties of GOE 
and lies at the basis of interpolation procedure widely used in the proof of Lemma 13.81 below. In 
fact this condition is pure technical one, and we can replace condition (11.21) with more general one 
and consider Wigner matrix M = n~ l l 2 W , satisfying 

E{^ } } = 0, 1 < j < k < n, (3.14) 

E{(W5 } ) 2 } = ™\ j ± k, B{(wjf) 2 } = w 2 w 2 , w 2 > 0. 

In this case there arise additional terms in (13.1ip and (13. 13ft proportional to w 2 — 2. In particular, 
we have for the corresponding limiting variance 



r2w 

V d w [cp] = V d w [tp] + (w 2 -2)w 



2w 



(3.15) 



where VW [cp] is given by (|3. 13|) . 



d 

Proof. Let us write the covariation in the form (cf (|3,8p 



nCov{(<p 1 (M)) jj ,(<p 2 (M)) jj }= I j C n {t u t 2 )$ x {tx)(p 2 {t 2 )dtxdt 2 , (3.16) 

where 

C n (ti,t 2 ) = nE{U jj (t 1 )U° j (t 2 )}, 
and U is defined in f|2 . 9[) . It follows from the Schwarz inequality and (|3.9p that 

CW(*i,*a) < C7(l + |ti|) 3 (l + |t 2 |) 3 - 
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This, (|3.4p . and f|3. 16[) imply that it suffices to show that there are converging subsequences {C n .} 
and function Cov such that we have for any converging subsequence {C ni } and any T > 

lim CntihM) = Cov(h,t 2 ) (3.17) 

i— >oo 

uniformly on the square St = { (t 1 , ^2 ) £ ^ : |^i| < T, \t 2 \ < T}, and that plugging Cov in the r.h.s. 
of ([3T61) we get the r.h.s. of (jSTTTI) . 

Show first that the derivatives dC n /dti, i = 1, 2, are bounded on St uniformly in n. We have 

^ = inE{(ME0;;(ti)E£-(t 2 )} = ^^E^^^i,^)}, (3.18) 

k=l 

where 

*ifc(ti,t2) = ^ifc(«i)t^(t 2 ). (3.19) 
A simple algebra based on (|2.15p allows to obtain 

\D l jk <f> jk {t 3 ,t 2 )\ = 0(l), n^oo, (3.20) 

uniformly in (ii , £2) £ St- Now applying differentiation formula (|2.20p with = <!>,,■£, and p = 2 to 
every term of the r.h.s. of (|3,18p . we get 

^ = iw 2 J2^E{D jk ^ jk (h,t 2 )} + ^J2 E i D %^(ti,t 2 )} + e 2 {t 1 ,t 2 ), (3.21) 

1 k=l ^ U k=l 

where in view of (HOT]) and (ET201 

\e 2 (h,t 2 )\ < V sup \D%$ jk (t lt t2)\ = 0(1), n-»oo. (3.22) 

Here 5 n is the set ofnxn real symmetric matrices. Now it follows from (|2,12p and (I3.19P that 
every term of Dj k &j k contains Uj k (see (|3.35p ). Taking into account that X]fc=i Wjk\ < n 1 ^ 2 ( see 
(12. lip ), we see that the second term on the r.h.s. of (I3.2ip is of the order O(l), n — > 00 uniformly 
in (ti,t 2 ) £ St- At last, using (|2.12p we get for the first term on the r.h.s. of (I3.2ip : 



k=l 
n 

■ 2 Y, E{(^ii * Ukk + U jk * UjkXh^ih) + 2U ]k {t 1 ){U j] * U jk ){t 2 )} 

k=l 

E{ [{Ujj * nv n )(h) + hUjjih)] U°At 2 ) + 2 I 2 U^h + U)Ujj{t 2 - U)dt 4 }, 

Jo 



iw 2 ' 1 



where we denote 

n 

v n {t):=n- 1 TxU{t) = n- 1 Y J U k k{t), \v n (t)\ < 1. (3.23) 

k=l 

Since 

nEivnUjjU^} = B{v n }C n + nE{<?7 ii £7?.} 
and by (13.1061) - (13.1071) of Lemma E31 below 

Var{^} = 0(n _1 ), Var{v n } = 0(n~ 2 ), n 00, (3.24) 
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uniformly in (t\,t 2 ) £ St, then we finally have 

Tl n \h,t 2 ) =iw 2 [ 1 E{u n (t 4 )}C n (ti - t 4 ,t 2 )du 

Jo 

+ 2iw 2 / EiUjjih + U)}E{Ujj(t2 - U)}dU + r n (ti,t 2 ), (3.25) 

where 

r n (h,t 2 ) =iw 2 E{ [(Ujj * m£)(ti) + C^(t 2 )} 

+ 2 / 2 E{?7 7i (ti+t4)^° ; (t2-t4)}^4 = 0(n~ 1/2 ), n-fr-oo, (3.26) 
jo 

and see that T[ n \t ll t 2 ) = 0(1), n — > 00 uniformly in (ti,t 2 ) £E <Sr- 

It follows from the above that the derivatives dC n /dti, i = 1, 2, are bounded on St uniformly in 
n. Hence, there are converging subsequences {C n .} and function Cov (depending on subsequence) 
such that (|3,17p holds. Now we derive an integral equation for Cov showing that Cov is the same 
for every converging subsequences {C n .} and leading via f|3. 16[) to p. lip . 

It follows from the Duhamel formula (12.81) that 



C n (h,t 2 )= I m 1 / 2 ^E{^ ) $ ifc (t 3 ,t 2 )}^3 



k=l 



with &jk(t3,t 2 ) given by ()3. 19j) . We see that the integrand here coincides with the r.h.s. of (I3.18p . 
Hence, applying differentiation formula (|2.20p with p = 3, we get (cf (|3.2ip - (|3.22p ): 



dt 3 , (3.27) 



C n (h,t 2 )= I' i y2Tl n \t 3 ,t 2 ) + e 3 (t 3 ,t 2 ) 
Jo I l=1 J 

where 

1 n 

T J (*3,*2) = Un (i-i)/2 E^l^^fe.^)}, I = 1,2,3, (3.28) 

fc=i 

(n) 

Kljk is the Ith. cumulant of Wj k : 

Kijjfc = 0, K 2j fc = w 2 /^ 1 , K 3,ife = M3, K 4)jfc = k 4 - ^a^ 4 , (3.29) 

(see CCD, (|2~T8D . and (13321) ). and in view of (123TD and (1330]) 

5/6 n 

k 3 (*3,*2)| < -^|- y; sup |D4$. fc ( t3 ,t 3 )| =0 („-i/2) > (3.30) 

We see that T 1 {n) of (133HD is given by ([3351) - (13351) . where by (I3.105P - (|3.1U6I) of Lemma 1331 
below 

lim E{77 n (t)} = lim E{[/ 7J (t)} = u(t) = / e itx p sc (X)dX. (3.31) 

n— >oo n— >oo J 

This and (l3~T71) yield 

lim r 1 (ni) (t 3 ,t 2 ) = iw 2 [ 3 v(t 4 )Cov(t 3 - h,t 2 )dt 4 + 2iw 2 $(t 3 ,t 2 ), (3.32) 
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where 



2w p2w g^2^ g^2M 



Ht 3 ,t 2 )=i- 1 I / e it3X p 8C (X)p sc (ti)dXdfi. (3.33) 

-2w J-2w A — fJ, 



Consider now T 2 (n) of (l3~28l) . and show that 



lim T 2 {n) (t 3 ,t 2 ) = 0. (3.34) 



We have by (f2~T2l) and (l3~28l) with 1 = 2; 



T 2 {n) (t 3 ,t 2 ) = —^LJ2P%V{(U jk * U jk * U jh + 3E^- * U jk * U kk ){t 3 )U°j(t 2 ) 



+ 2([/ jfc * [/ jfc + Ujj * U kk )(t 3 )(Ujj * U jk )(t 2 ) 

+ U jk (t 3 )(3Ujj * U jk * U jk + Ujj * Ujj * U kk )(t 2 )}. (3.35) 

It follows from (|2.1ip that the contribution of the terms containing Uj k Uj k Uj k is of the order 
0(n -1 / 2 ), n — > oo. Besides, since n~ l l 2 Y^k=i \U k k(t)Uj k (T)\ < 1, then by the Schwarz inequality 
and (|3.9p the contribution of the terms containing Ujj(t 2 ) is also of the order 0(n~ 1//2 ), n — > oo. 
So we are left with 



— % E E i 2 (^ * Ukk)(t 3 )(Ujj * U jk )(h) + UjkfaXUjj * Ujj * U kk )(t 2 )}. 

n k=l 

Here by (13.1091) of Lemma EH] below 



n 

Varjn- 1 / 2 ^ U jk (n)U kk (r 2 )} = O^- 1 / 2 ), n -> oo, (3.36) 
fc=i 

lim Ejn- 1 / 2 V C/ ife (ri)C/ fefe (r 2 )j = 0, (3.37) 



fc=l 



so that 



n 

lim Ejn- 1 / 2 V ^(71)^(74)^(73)^^4)} =0, 



fc=i 



and we get (|3.34p . 



Consider now T 3 (n) of (pi28]) . We have 



n 



T 3 (n) (t 3 ,t 2 ) = -^^E{D| fc (C/ ife (i 3 )^(i 2 ))} + 0(n- 1 ; 

n k=l 

where we replaced K4j k of (|3.29p with K4 of (|3.12p with the error term of the order 0(n -1 ), n — > 00. 

r (r ' 



It follows now from (|2.1ip - f|2. 12j) that the contribution to Tg due to any term of 



n 

— 1 ' 

n 

k=i 



^D^Uj^U^h)) 



containing at least one off-diagonal element Uj k is of the order 0(n 1 ^ 2 ), n — > 00. Besides, we have 
by (I2J51) and ([32]) that the term E{U°j{t 2 )n- 1 Y2=i D 1k U jk{^)} is of the order C^n" 1 / 2 ), n -> 00, 
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too. Thus, we are left with terms, containing only diagonal non-centered elements of U. There is 
only one such term, it arises in the term 3DjkUjk(t 3 )Dj k Ujj(t 2 ) of the sum above, and by (I2.12P its 



(n) 

contribution to T 3 is 



~ V Yl E U u 3i * u kk){h){u 3J * Ujj * u kk )(t 2 )}. 



k=l 

In view of (|3.9p . ([3.106p . and (|3.108p we can replace here all Ujj and U k k with v in the limit n — > oo, 
so that we have 

lim To (n) (t 3 ,t 2 ) = -K4i(v*v)(t 3 )(v*v*v)(t 2 ). (3.38) 

n— >oo 

Summarizing (|3,27p . (|3.30p . (|3,32p - f|3.34|) . and (|3.38p we obtain the equation with respect to Cov 
of (13471) : 

rh rta 

Cov(t u t 2 ) + w 2 dt 3 v(t A )Cov(t 3 -t 4 ,t 2 )dU = A(h,t 2 ), (3.39) 
Jo Jo 

where 

A{h,t 2 ) = -2w 2 / ${t 3 ,t 2 )dt 3 + k 4 (v * v * v)(t 2 ) / (v*v)(t 3 )dt 3 (3.40) 



and <E> is given by (|3.33|) . To solve (|3.39p we use the generalized Fourier transform with respect to 
t\ (see Proposition 12. ip . Note, that equation f|3.39[) is of the form f|2.4[> . corresponding to 5 = in 
(12.11). thus we can use formulas (12.51) - (12.61) to write its solution. Since 



v(z) := -i / e- itz v(t)dt = {2w 2 )- 1 {^Jz 2 -Aw 2 - z), (3.41) 



with the branch that is determined by the asymptotic \J z 2 — 4w 2 = z + 0(z 1 ), z — > oo, and 
(z + w 2 v(z))~ 1 = —v(z), then we have for T of f|2 . 6[) : 

T(t) = — I e ltz %— - = -v(t). (3.42) 

V ' 2vr J L z + w 2 v(z) W V 1 

Hence, the unique differentiable solution of (|3.39p is given by 

rti 2 
Cov(h,t 2 ) = -2w 2 / v(ti -t 3 )<S>(t 3 ,t 2 )dt 3 + K 4 Y[(v*v*v)(tj). (3.43) 

jQ 3=1 
We have by (pT33|) . (|3.105l) . and a little algebra: 

/ V(t 1 -t 3 )$(t 3 ,t 2 )dt 3 = - T[p sc (Xj)dXj, 

JO J\-2w,2w]' i M — M M — A3 fj~ 



and by (f3T4"Tj) 



(«*«*<;)(*) = — / e ^(2u» 2 )- 3 (^z 2 - 4u> 2 - z) 3 cfc 

27T 



' L 

r2w 

w~ 4 I e itx (w 2 - X 2 )p sc (X)dX. 

-2w 
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Hence, putting these expressions in f|3.43|) . and then plugging the result in (|3,16p . we finally get: 



Jkn nCov{ (ipx (M) ) jj , (<p 2 (M) ) jj } 



2 [ (V>i A i ~ fi(M))(<P2 A 2 ) - y 2 (A 3 )) t r / \ x , x 

J[-2«),2«)]3 (.Ai - A 2 J(A 2 - A 3 J ^ 



2 /-2ji! 

+ K4 J] / ^(A)(^ 2 - A 2 )p sc (A)dA. (3.44) 



<"2mi 

Writing the numerator in the first integral as 

^i(Ai)y?2(A 2 ) - ¥>i(Ai)<£ 2 (A 3 ) - 9?i(A 2 V2(A 2 ) + <Pi(A 2 )v 3 2(A 3 ) ! (3.45) 
we observe that there is at least one integration which does not involve (p's. This and the relation 

p sc (p)dp x /0 __ 2 



p — X 



allow us to deduce (|3.1ip from (I3.44p . A simple way to perform the corresponding calculations is 
to write the r.h.s. of (|3.44p as the limit as e — > of the same expression in which A 3 is replaced 
by A 3 + is. One can also use the Poincare - Bertrand formula [10J to deal with double singular- 
integrals, appearing after plugging (I3.45P in t!3.44j> . ■ 

3.1 Limit Theorem for Matrix Elements 

Theorem 3.4 Consider the real symmetric Wigner random matrix of the form 

M n = n-^Wn, W n = {W jk € R, W jk = W kj = (1 + 5 ifc ) 1/2 ^ fe }£ fe=1 , (3.46) 

where {Vj k }i<j< k<OQ are i.i.d. random variables such that 

B{V n } = 0, E{Vfi} = w 2 , 

and the functions lnE{e li1/l1 } are entire. 

Then for any (p : R — > M whose Fourier transform \3. 3\) satisfies tfTM and for any j = j n E [i-,n] 
the random variable y/nip^j^M) converges in distribution as n —> oo to the random variable £ 
having the characteristic function 

E{e ix S} = exp{{-x 2 V^[ip] + w 2 x* 2 )/2}f(x*), (3.47) 

where f(x) = E{e lxVl1 }, 

x* = — T I (p(p)pp sc (p)dp, (3.48) 

w J-2w 

p sc is the density of the semicircle law hi. 7|) , and V^r [tp] is given by tS.13)) . 

Remark 3.5 Condition W jk = (1 + 5 jk ) l/2 V jk is pure technical. In particular, it can be shown 
that in the case of matrix M = ra~ 1//2 V ', the Theorem 13.41 holds true with 



w 2 



x f 2w 

2w 



and 

E{e ia *} = exp{(-a^V^[ip] + 2w 2 x* 2 )/2}f{a 
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Proof. Note first that in view of (|2.14p and (|3,46p we can write 

W jk = PJk /2 Vjk- (3.49) 
Besides, since lnEje**^ 11 } is entire then we have 

Y,^JT^<oo, Vx>0, (3.50) 
1=1 L 

where k\ is the Zth cumulant of V\\. We also have 

ln/(x) = -^V/2 + £^M- (3.51) 

1=3 L 

(see pTfj) ). 

We consider the characteristic functions 

Z jn {x) = E | e -^(M)| (3 52) 

and prove that for any i£l 

lim Z jn (x) = E{e ia *} =: Z d {x), (3.53) 

n— >oo 

i.e., Zd(x) is given by the r.h.s. of (|3.47p . 

Assume first that the Fourier transform fj3.3j) of satisfies 



|£(i)||i| z ctt < Cy! VIeN, (3.54) 

in particular, ip is analytic in \z\ < 1. Since Zj n (0) = 1 and Zj n (x) is continuous, we can write the 
relation 

Z ift (x) = 1 + jT* Zj n (y)dy, x£R, (3.55) 

showing that it suffices to prove that the sequence {Zj n } is uniformly bounded on any finite interval 
and that for any converging subsequences {Zj ni }i>\ and {Zj n .}j>i there exists Z(x), such that 

lim Z jni (x) = Z(x), (3.56) 

i— >oo 

and 

Hm Z' jn .{x) = Z{x) - xVf\ip\ + 7fZ TVj (~^T [ pM^cM^) • (3.57) 
Indeed, if yes, then Z{x) is a continuous function, satisfying for every i£l the equation 

Z{x) = l- j* Z{y)[-yVfM 

J-l , iy /2 r2 W 



(7_ i)! (^T / 2 vQ^WscQ^dlt) dy, (3.58 



whose unique solution is the r.h.s. of (|3.47p . 

We denote _ 

^ n (x) = e* e ^W, (3.59) 
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and write according to (|3,7p and (13.5211 

Z' jn (x) = iE [^ 3] {M)e ix ^°n {M) ] = i J <p(t)Y jn (x, t)dt, (3.60) 

where 

Y jn (x,t) = v^E{^(t)e° n (x)}, (3.61) 

and U is denned in (I2.9p . It follows from the Schwarz inequality and (|3.9p that 

\Y jn (x,t)\<C(l + \t\) 3 . (3.62) 

This and (I3.54p with I = 2 yield that the sequence Z'- n is uniformly bounded. Hence, there is a 
convergent subsequence Zj n i , and by the dominated convergence theorem to find its limit as n — > oo 
it suffices to find the pointwise limit of the corresponding subsequence Yj n / . 

Let us show now that sequences {dYj n /dx} and {dYj n /dt} are uniformly bounded in (t, x) G 
K C M+, n G N, for any bounded K, so that the sequence {Ijn} is equicontinuous on any finite set 
of and contains convergent subsequences. 

Since Yj n (x,t) = Yj n (—x, —t), we can confine ourselves to the half-plane = {t > 0, x G M}, 
and from now on t > 0. 

It follows from (|3.7p that 

—Y jn (x,t) = i / ^(t 1 )nE{C/° ? .(t 1 )C/? ? .(i)e iri (x)}^i, 

where by (|3.9p and the Schwarz inequality 

n|E{C^.(ti)E/?,.(t)e in (x)}| < nVar 1 /2{i/ ii (t 1 )}Var 1 /2{ E / ii (t)} < C(l + |t|) 3 (l + |^|) 3 . 

Hence, in view of (|3.54p the sequence {dYj n /dx} is uniformly bounded. 
We have also 

a n 

g-Y in (x,t) = i^V{{MU) 3j {t)e° jn (x)} =iJ2 E {W jk $jk(x,t)}, (3.63) 

k=i 

where 

<f> jk (x,t) = U jk (t)e° n (x). (3.64) 

To transform the r.h.s. of (|3,63p and show its boundedness, we apply an analog of integration by 
parts formula proposed in Lemma 12.31 Note that Dj k &jk = 0(n^ 2 ) as n — > oo, hence, there is no 
such finite p G N that e p of (I2.20P vanishes as n — > oo, and so we need infinite version of "integration 
by parts formula" given by (|2,22p . We will apply f|2.22[) to every term of the r.h.s. of (|3.63p . and 
to do this we check first that &j k (x,t) satisfies condition (I2.23p . Indeed, using the Leibnitz rule we 
obtain 

I , 

D l jk $ jk (x,t) = J2 ( m )D l jk m U jk (t)D'? k e° n (x), (3.65) 

m=0 

where 

Df k e jn {x) = DJ- 1 (ix^e jn (x)D jkVjj {M)) , (3.66) 

(see (j3~59~j) ), so that 

m r 

Df k e 3n {x) = e jn (x) ^(ix^) r £ C q , r ]] D% Vjj (M), 

r=1 q=( qi ,...,q r ): s=1 

qi + ... + q r = m 
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and 

E Cq,r < 2 m . 
q,r 

Hence, 

r 

\Df k e jn (x)\<{2^l(l + \x\)) m max \{\D%^{M% 

where we have in view of f|2. 15[) and ()3.54[) 

\D%^{M)\ < J \m\\D%U n (9)\de < C V V>, (3.67) 

so that 

\D™ k e jn (x)\ < (4C^(1 + \x\)) m - (3.68) 

This, and (^65]) yield 

|D} fc * ifc (x,t)| < (4C7^>/n(l + N xGM, t>0. (3.69) 

Thus, $>jk(x,t) satisfies condition (I2.23|) . Applying (12.22P to every term of the r.h.s. of (I3.63P and 
taking into account (|3.49p . we get: 

|y iB (x, t) = < E ' s ™ ( *' t} = E ^^^{^(x, *)}■ (3-70) 

i=i ' w ^ fe=i 

Let us show that this series converges uniformly in (t, x) £ K, n G N. In view of ([3.50P it suffices 
to show that 

\s\ n) \ < (C K ) 1 , V(t,x) €K, n G N, (3.71) 
where Ck is an absolute constant depending only on K. Since 



2 (l+l)/2 _ i 1 ™_ 



fc=i 



(see (|2.14p ). where in view of (|3.69p the first term of the r.h.s. is bounded, it suffices to prove (|3,7ip 
for the second term of the r.h.s. of (I3.72p . 
Using the Leibnitz rule, we write for I > 2 

1 n 1 n I n 

J^y E D j^Mx, t) = ^ E U ^D l jk e jn + ^ E D jk U jk D l ^e jn (3.73) 

k — 1 k — 1 k — 1 

1 n 1-2 

+ T-W E E ( m )^U, k D^ n =: atf + a&> + «#>, 

VV ' fc=lm=0 

where (cf (|3~65l) - (j!TB5|> ) 



|a[ 3 n) | < (4(7^(1 +>|+t))'. (3-74) 



Applying f|3.66H and then the Leibnitz rule again, we obtain for a,^ of f|3.T3[) : 



(n) «c 



^ E U ]k {t)D jWn {M)D l r k l e jn (3.75) 



n Z-2 . _ 

vv 7 fe=l m=0 
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where the sum over m is bounded by 



(y/E) l - 2 {4C v (l + \x\+t)) 



l-l 



(cf (|3.69j) and (j3.74D ). Taking into account that \ "Yuk=iUjk(f)\ — n 1 ^ 2 > we see that the second 
term of the r.h.s. of (|3,75p is bounded by 

(8CV(l + |x| +*))'. 
Besides, it follows from ([221) and flZ3Z]) that 

D jk<Pjj (M) = 2i/3 jk f tp(0) J° Ujjip - e^Ujk^dBidO, (3.76) 

roc rd 

\D jWjj (M)\ < 2 / / |*7 ifc (0i)|d0id0 < 2C V . (3.77) 



This, (|3.67p - (|3.68p . and (|2.1ip allow us to show that the first term of the r.h.s. of ([3.75P is bounded 
by 



2\x\(4C v (l + \x\)) 1 1 I \tp(9)\ I Y,\ U #®W u M i)\Mi M 



k=l 



< (4c v (i + \x\ + \t\)y. 



Hence, 



An) 
hi 



< (4C v (l + \x\+t)) 1 . 



(3.78) 



Finally, since by (l2~T2l - (l2~14l 

Dj k U jk (t) = i(Ujj * U kk + U jk * U jk )(t) - i5 jk (Ujj * Ujj)(t), 

we have for a™ of (13.73|) : 



>) 

hi 



a 



u 



, £(^i * U kk ){t)D l - k l e ]n + T7 ^- l J2(U jk * U jk )(t)D l r k i e 

W ' k=l ^ V ' k=l 



il 



j^(U n *U n )(t)D^ ejn 



where the last two terms are bounded by 1(1 + i) (4(7^,(1 + |x|))' 1 in view of (I3.68p . (12. lip , and the 
bound \(U a b*U c d)(t)\ < \t\. Besides, it follows from (|3.66j) and (|3.77p that the first term is bounded 
by: 



It 



Efe 1 ' 



1-3 



k=l 



" T Y, \D jkiPjj D^e jn • E ( !., 2 )D l f k 2 - m ^D Tk e JT 



< 



k=l 
lt\x\ r 



m=0 



in 



-i 



n 



2 Y / \u jk (e 1 )\de 1 d0(AC (p y/E(i + \x\)) 

Jo Jo k=1 



1-2 



+ (V^y- 1 (4c v (i + \x\ + \t\j) 



(1-2) 



< (4C v (l + \x\ + t)) 1 



so that 



l l2 I 



< (40^(1 + \x\+t)) 1 



(3.79) 
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Now (HTTP with / > 2 follows from (pTT2]) - ([3~T4l and ([3~T8D - (pT79]) . Hence, the series in (l3~63l) 
converges uniformly in (t, x) £ K C and n £ N. 

To prove the boundedness of the sequence {dYj n /dt}, it remains to make sure that is 
bounded. Applying fl2~T2j) - (l2~14l) and fETKHI) . we obtain 

n 

S[ n) (x,t) = n- 1 / 2 J2^MD jk ^k(x,t)} 
fe=l 

= m-^Ej (n(Z7^- * u n )(i) + tU j:j {t))e jn {x)} 
- 2x J 0(9) &{U 3j (9 - 9 x )U 3j (t + Qx)e jn {x)}de x dQ (3.80) 

where v n is defined in (|3.23p . Writing 

n 1 / 2 E{C/ ij (t 1 )7; n (t 2 )e° n (x)} 

= V{v n {t 2 )}Y jn {xM) + n 1/2 E{C/ i ,(t 1 )t>°(t 2 )e° n (x)}, (3.81) 

and taking into account bounds f|3.62|) . (|3.107p . \Ujj\ < 1, and |e? n (x)| < 2, we conclude that the 

r.h.s. of (|3.81|) is bounded, and so does . Hence, the sequence {dYj n /dt} is uniformly bounded 
in (t,x) S K C R\, n E N. 

Now it follows from the above that the sequence {i^n} is equicontinuous on any bounded set of 
M 2 . Hence, for any converging subsequence {Zj ni } (see (|3.56p ) there is a converging subsequence 
{Yj n /} and function Y (which obviously depends on {Zj ni }) such that 

lim Y jn , = Y, lim Z jn , = Z. (3.82) 

We will show now that Y satisfies certain integral equation leading through (|3.60p to (|3,58p . hence, 
to (|3,47p . This will finish the proof of the theorem under condition (|3.54p . 
Applying the Duhamel formula (|2.8p and then (|3.63p and (|3.70p . we obtain 

ft n ft 00 

Y jn (x,t)=i^E Y j K{M jk <S> jk {x,t l )}dti = i V^sf n) (x,ti)dti, (3.83) 
Jo k=i Jo 1=1 11 

where and Sj are defined in (|3.64p and (|3.70p . respectively. In view of the uniform convergence 
of the series, to make the limiting transition as n — > oo it suffices to find the limits 

Si = lim 5, (n) 

for every fixed I G 1. 

Let us start with s[ n) . It follows from flcTHO]) - (l3~8TH that (cf (l3~25l) - fl3~26l) ) 

5| n) (x,ti) =i / 1 B{v n (t 1 -t 2 )}Y jn (x,t 2 )dt 2 (3.84) 
Jo 

- 2xZ jn (x) J 0(9) J V{U 33 (9 - 6 1 )}E{U 33 (t 1 + BifidBxdB + r n (x,h), 

where 

r n (x, t) =m- 1 / 2 E{ [n(Ujj * v° n )(t) + tU 3j {t)] e° jn (x)} 

-2x\ 0(9) [ E{U 33 (9 - ${) [U°j(t + 9 1 )Z jn (x) + U j3 {t + 9 1 )e° jn (x)] }d9 l d6, 
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and by (I3.24p and boundedness of U a b and ej n we have 

r n = 0(n-^ 2 ), n-+oo. (3.85) 

This, (13T3T|) . and (^82|) imply 

lim S[ n ' l) =i [ 1 v(h - t 2 )Y(x,t 1 )dt 2 - 2xZ{x) [ !p(d)^(t 1 ,6)de, (3.86) 

where $ is given by (|3.33p . 

In general case we have for sf n) , / > 2, of (I37DD : 

1 n 

S \ n) = T7^\i E l 3 jk +1)/2E { u jk Dl jk e jn + lD jk U jk D l ~ l ej n 



k=l 

1-3 



+ 1(1 - l)D 2 k U jk D l r 2 e jn /2 + (1 - 6 l2 ) £ ( ^ )D l r k m U jk Df k e jn \ 

m=0 



=4 n) + 4° + 4 n) + C 1 - w ( ? > ( 3 - 87 ) 

and we have (cf (13.74P ) 

\S% ] \ < (4^(1 + |x| + |i|)) V 1 / 2 , / > 3. (3.88) 

Now we use the rule that in fact has been used several times before and which follows from 
(12. lip and the boundedness of ej n of (13.59P : the presence of a single factor U ak in terms of the sum 
Y^k=i i s equivalent to the presence of the factor n -1 / 2 , and the presence of two or more factors 
U ak , Ub k ,... is equivalent to the presence of the factor n~ l . It follows from (I2.12P that all terms of 
Dj k Uj k contain Uj k , besides, we have 

D jk 2e jn = 0((V^) l ~ 2 ), n^oo (3.89) 

for I > 2, hence, 

^3 

We also have (see (I3.75P ): 



5 z (n) = 0{n~ 1 / 2 ), n^oo. (3.90) 



+ (/ - l)Dj k ipjj (M)D l j k ~ 2 ej n ] } + Oin- 1 / 2 ), n^oc, 



k=l 

+ l)D 2 k tpjj(M)L- k 
where Dj k (fjj(M) is given by (|3.76p . and 

D% Vjj {M) = -2i(3 2 k J ${e){W j3 * U jk * U jk + U j:j * U 3j * U kk )(6)d6, (3.91) 
Hence, using again the above rule, we get 

/ n 
m ~Q^ E^fc E {M*i)(^ * U jk ){9)D)r k \ jn 

+ i(l - VUjkfaXUjj * Ujj * U kk )(6)D l - 2 e jn }d9 

+ 0(n" l/2 ), n^oo. (3.92) 
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It follows from ()3.36j) - (|3.37p that if I = 2 then the contribution of the second term of the r.h.s. of 
(|3.92p vanishes as n — > oo. To show this in the case I > 3 we use the evident relation 

D%e jn {x) = {ix^D lWjj (M)) m e jn {x) + 0({^i) m - 1 ), n ^ oo, (3.93) 

with Djk<Pjj(M) given by (13.76j) . and conclude that the second term contains either UjkUjk ■ 
0{(y/n) l ~ 2 ) or Ujk ■ 0((^/n) l ~ 3 ), hence, we have in view of (|2.1ip : 



n -(J-l)/2 fikUjkihXUjj * Ujj * U kk )(6)D l r k 2 e jn (x) = C^n" 1 / 2 ), n ^ oo. 

k=l 

The first term of the r.h.s. of (13. 92ft already contains Uj k Ujk- Thus its non-vanishing is due to the 
term {ix^/nDj k ipjj{M)) l ~ 1 ej n {x) of D l ^}ej n {x), and we get: 

n 
k=l 

= Y,E{e jn {x)U jk (t 1 )(-2x / ip(0)(U jk *U j:j )(6)d6\ } 
k=i 

+ ( 2 (i-0/a _ l)E{e in (x)C/ ii (t 1 )(- 2x J tp(9)(U 3j * U j3 )(6)de) 1 } 
+ 0(riT 1/2 ), n-too, 

where we took into account (f3776|) and the equality ^ 1)/2 = 2( 1_ 0/ 2 ( see (|2TR]) ). Applying (|3.104p , 
(I3.106P . and (13.1101) we get for / > 2: 



7i t 

k=l m=l 

i 

= Z jn (z)v n2 (ti,0i, ]J E {^i(^n)} + 0(n- 1/4 ) 

m=l 

i 

= Z(x)v(ti) f (6» m )u(r m ) + o(l), m -)• oo, 



m=l 



and 



E{e jn (x) [] L^(r m )} = Z(x) JJ u(r m ) + o(l), nj — > oo. 

m=l m=l 

Besides, we have in view of (|3.105p 

(v*v)(6) = ^ / j^wM*!*. (3-94) 

w J-2w 



The above allows to write 

iy/2r r2,r 

'- x : ' v w 2 J_ 2w 



lim = V2Z(x) V (t 1 )(^ / 2 " <p(riWsc(»)dn) 1 . (3.95) 
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It remains to analyze sffl of ([3.87p : 
/ n 

4° =J^E^t 1)/2 ^D jk U jk D^ ejn } 



k=l 

n 



Ix 



lX — E ^ 1V2e {(^' * U kk + U jk * U jk ){t{)D l r k 2 {e jn {x)D jWjj {M))} 



k=l 



71 



k=l 



+ (1 - <fe)(Z - 2) J D^ 3 e, n (x)Z)2 fc ^ i (M)) ) + O^" 1 / 2 ), n oo, 



where we used consequently (|2.12|) . (|3.66|) . and then the Leibnitz rule, (|3.68|) . and (I2.1ip . Treating 
the first term of the last expression analogously to the second term of (|3,92p . we see that it is of 
the order 0(n -1 / 2 ), n — > oo. Hence, taking into account (I3.91j) and (12. lip , we get 5*22 = 0, and for 
I > 3 

r i n 

S% ] = 2l(l-2)x J QWjj^Y.Efa^UkMh) 

x (Ujj * Ujj * U kk ){e)D l r* ejn {x)}dO + O^" 1 / 2 ), 



n — > oo. 



If I > 3, then in view of (T3793D . D l - 3 ejn (x) gives either factor U jk O((^/n) 1 - 3 ) or 0{(^i) 1 - 4 ), that 
in both cases leads to 



'12 

If I = 3, then 



Sjjp = 0(n -1 / 2 ), n^oo, I > 3. (3.96) 



5g } =6x / £(0)~^E{(% * UkkXhXUjj * ^ * [7 u )(0)e in (x)}d0 



fe=l 

+ <3(n~ 1/2 ), n-^oo, 

and it follows from (|3.1()2j) . (I3.106p . and (I3.108p that 

n 

E{e jn Ujj(Ti)Uj j (T 2 )U j j(T3)rr 1 ^ U kk (Ti)U kk (T 5 )} 

k=l 



Z jn (x) 11 BiUjjiT^Bin-^UkUr^U^)} + Oin- 1 / 2 ) 

m=l k=l 

5 

Z(z) u(r m ) + o(l), ni ->■ oo. 



m=l 

Hence, 



lim = &cZ(a:)(t;*t;)(ti) / £(0)(w * v * u)(0)d0. 



This, (I3T871) - fl3~90]) . and (l3~95]l - (l3~96l) yield for Z > 2: 

iy/2x r2w 

W 2 J- 2w 

+ 5 3l 6xZ(x)(v*v)(t 1 ) I (p(9)(v*v*v)(0)d0. (3.97) 



lim S\ ni) =Z(x)v(t 1 )( t -^^- [ <p(n)np 8C ((i)dii) 



19 



Summarizing f|3.83j) . f|3.86|) . and (|3.97p we see that Y of (|3,82p satisfies the equation 

Y(x,t)+w 2 [ dti [ v(ti-t 2 )Y(x,t 2 )dt 2 (3.98) 
Jo Jo 

= ixZ(x) / tp(6)A(t,6)d6 + iV2B(x)Z(x) [ v(h)dti, 



where A is defined in f|3.40[) . and 

B W = zZnr C-^ir / vMwMdn) ■ (3.99) 

]Zn W J~ 2w 

The kernel of this equation coincides with that of ([3.39p . Hence, the argument leading to (|3.43p 
and based on using of generalized Fourier transform, yields 

Y(x,t) = ixZ(x) [ fi(8)Cov(t,6)d8 + iV2B(x)Z{x)(v*v)(t), 



where Cov is given by (|3.43p . Plugging this expression in (|3.60p and taking into account (|3.94p and 
the equality (see (|3.16|) ) 

J j tp{6)<p{t)Cov{t,0)dtd6 = V] v [ip\, 

we finally get for Z of (^52]) : 

iV2 ^ 



Z'{x) = Z(x) 



w- 



2w 



This yields (|3.56p - (|3.57p . thus proves the theorem under condition (|3.54p . 

The case of ip G E = {ip : J (l + \t\) 3 \tp(t)\dt < 00} can be obtained via a standard approximation 
procedure. Indeed, since the set D = {tp : f \(p (t)\\t\ l dt < C^ll, V/ G N} is big enough (in particular, 

2 

it contains functions e~ x P m {x), where P m (x) is a polynomial), then for any tp G E there exists a 
sequence {<pt} C D, such that 



2ir 



lim / \cp(\) - <p k (\)\d\ = 0. (3.100) 

fc^oo J_ 2w 

Denote for the moment the characteristic functions of (|3.52|) and (|3.53|) as Z n [<p] and Z[ip], to make 
explicit their dependence on (p. We have then for any <p G E 

\Z n [(p] - Z[ip]\ < \Z n [p\ - Z n [ip k }\ + \Z n [ip k ] - Z[tp k ]\ + \Z[(p k ] - Z[ip]\ 

■^T^+tW+T®. (3.101) 

The second term of the r.h.s. vanishes after the limit n — > 00 in view of the above proof, since 
ip k G D. For the first term we have from (|3.52p and the Schwarz inequality that 

\T^\ < |a?|(nVar{(^(M)) ii }) 1/2 ) ip k = <p-<Pk, 
and then Theorem 13.21 implies that 

hmsupir^l < \x\(V^k}) 1/2 - 

n— >oo 

Since Vj^ of A3. 131) is continuous with respect to the L 1 convergence, then in view of (I3.100P T^j} 
vanishes after the subsequent limits n — > 00, k — > 00. 

At last, we have by ([3.50P and the continuity of the r.h.s. of (|3.47p with respect to the L 1 
convergence, that the third term of (|3,10ip vanishes after the limit k —> 00. Thus, we have proved 
the Central Limit Theorem under condition (13.41). ■ 
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Corollary 3.6 It follows from Theorem \3.4\ that if ip is even, then the random variable y/n(f°(M)jj 
converges in distribution to the Gaussian random variable with zero mean and the variance VW [<p]. 

Remark 3.7 Random matrix theory deals mostly with eigenvalues of large random matrices. How- 
ever, the statistical properties of eigenvectors are also of considerable interests for a number of rea- 
sons, in particular in view of possible links with the problem of existence of absolutely continuous 
spectrum of the multidimensional Schrodinger operator with random potential (see e.g. (3j [5]). In 
the case of the Gaussian random matrices (GOE, null Wishart) the eigenvectors are rotationally 
invariant and according to recent works [U EJ [6] the eigenvectors of the non-Gaussian random ma- 
trices (Wigner, sample covariance) are similar in several aspects to the eigenvectors of the Gaussian 
random matrices. On the other hand, the results of [8] and this papers imply that there are asymp- 
totic properties of eigenvectors of the non-Gaussian random matrices which are different of those 
for the Gaussian random matrices. 

3.2 Auxiliary results 

Lemma 3.8 Consider the unitary matrix U(t) = e ltM of K2. 9\) - i\2.10\) . where M is the Wigner 
matrix M-l\) - M.2\) . define 

n 

Vfl (t) =n~ 1 ^2ll kk (t), 

k=l 

n 

v n (ti,t 2 ) =n- 1 Y,U k k{ti)U kk {t 2 ), (3.102) 



fc=i 

n 

«nl(*l,*2) =n- 1/2 Y, U Mti)U kk {t 2 ), (3.103) 
k=l 

n I 

Vn2(t) = II U lk{tm), I > 3, (3.104) 

k=l m=l 



where t = (t\, -.,ti), and put f = E{/}. Then we have: 

e itx p sc (X)d\ =: v(t), (3.105) 

-2w 

where p sc is the density of the semicircle law (il. 7\) . and under the conditions of Theorem \3.2\ 

(i) V&riUjjit)} < C(l + \t\f/n, lim U j:j (t) = v{t), (3.106) 

(ii) Var{v n (t)} = 0(n" 2 ), n oo, (3.107) 

(iii) Var{v n (ti,t 2 )} = 0(n- 1/2 ), n ^ oo, lim v n {t x ,t 2 ) = v(t!)v(t 2 ), (3.108) 

n— >oo 

(iv) Var{v nl (ti,t 2 )} = 0(n" 1/2 ), n -> oo, lim u„i(ti, t 2 ) = 0, (3.109) 

n— >oo 

I 

(v) Var{v n2 (t)} = Oin- 1 / 2 ), n ^ oo, lim v n2 {t) = Y[ v(t m ), (3.110) 

m=l 

where 0(n a ) can depend on t, t. 

Remark 3.9 It can be shown that all statements of the lemma remain valid under conditions of 
Theorem 13.41 
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Proof. Statement (|3.105p follows from the well known fact of random matrix theory (see e.g. 
12] and references therein) according to which for any bounded and continuous ip 



lim n _1 E{Tr cp(M)} 



2w 



ip(X)p sc (X)d\, 



2w 



where p sc is the density of the semicircle law (jl.7D . 

(i) Let M = rT l l 2 W be GOE matrix JO} independent of M, and 



itM 



We can write 



U(t) = e 



V n := VariUjjit)} = B{(U n (t) - U^tyU^-t)}, 



(3.111) 



(3.112) 



and then follow the interpolation procedure proposed in [8]. Namely consider the "interpolating" 
random matrix (see |13^ [8]) 



M{s) = s 1/2 M + (1 - s) 1/2 M, 0<s<l, 



(3.113) 



viewed as defined on the product of the probability spaces of matrices W and W. We denote again 
by E{. . . } the corresponding expectation in the product space. Since M(l) = M, M(0) = M, then 
putting 



U(t,s) 



JtM{s) 



we obtain 



u jj (t)-u jj (t) = y o fu n (t, s )ds 







(3.114) 



(3.115) 



l n 



Lm=l 



w, 



(n) 



Irn 



V(l - s)n 



Wi m )(Uji * U mj )(t,s)ds. 



Thus, 



■ fl r i n 1 n 

Lm=l v V > l,m=l 



ds, 



where 



3>lm = (U j l*U mj )(t,s)U° j (-t). 
A simple algebra based on (I2.12P and (12.15^ allows to obtain 

\D? m <s> lm \ < c q \tr\ 



(3.116) 



(3.117) 



with C q depending only on q £ N. 

Now, applying differentiation formula (|2.20p with p = 4 and $ = <^ m to every term of the 
first sum and differentiation formula (|2,16p to every term of the second sum, we obtain (cf (I3.27P - 



Vn 



+ £4 



- x / 2 ds, 



(3.118) 
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where 

n 



T v n) = -^m E ^E{C^}> P = 2,3,4, (3.119) 

i,m=l 

and by (HOB and (I3.117P 



M < V sup |Z>f m * Im | < C(l + It^n- 1 . (3.120) 

n i, m =i MeS " 

Consider now T 2 n ^ and note that by (|2. 18|) and (|3.ip K3 j m = 1U3, so that 

n 

T ^ = E E{A 2 AI (3.121) 

Z,m=l 
n 

= 2^ E V{U° 3J {-t)D* m (U 3 i*U m3 ){t,s) 

l,m=l 

+ (U jl *U mj )(t,s)D? m U JJ (-t) 

+ 2D lm (Uji * U mj )(t, s)D lm U jl (-t)} := M3 [T 2 ( ? + T 2 ( 2 n) + T^]. 

It follows from ([2TT2]) that 

n 

T 2 ( x n) = -3^ 3 S^ 3/2 E E {^(-*)(^ * * * C/Zm (3.122) 

i,m=l 

+ 2Ujl * Ujl * Ul m * Umm + Ujl * Uj m * Uft * U mm )(t, s)} . 

Here by (|2~TTD 

n n 
E PilUjmVM < 1, PjlUjlUimUmml^n 1 / 2 , (3.123) 

l,m=l l,m=l 
n 

\UjiU jm U u U mm \ < n. (3.124) 

l,m=l 

Hence, applying the Schwarz inequality and taking into account that (1 * 1 * 1 * l)(i) = i 3 /6, we 
obtain 

\T^\ < C(l + \t\ 3 )(n-V 2 Vy 2 + n- 1 ). (3.125) 
The terms T 22 ^ and r 2 g contain sums of two types 

n n 

UflUjiUjtUjmUnm = 0{n 1 ' 2 ), £ U^U^U^U^ = O(l), n ^ oo, (3.126) 

Z,m=l Z,m=l 

where the r.h.s. of both equalities follows from (|2.1ip . Hence, |T 22 + < C(l + |t| 3 )n _1 . This, 
([3TT2T]) . and (pm5T) yield 

|T 2 (n) | < C(l + |i| 3 )(n- 1/2 F n 1/2 + n- 1 ). (3.127) 

Acting in the similar way and taking into account (12.19P and (13 .21) implying < C, p = 4, 5, 

we get analogous bounds for T 3 (n) and T 4 (n) of (l3~TT8l) : 

|r/ n) | < C(l + |t|) ,+1 n _1 , Z = 4, 5. (3.128) 
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Note, that in the case of and the argument is even simpler because here we have factors 

n~ 2 and ra~ 5 / 2 , respectively, instead of n~ 3 / 2 of T^. 

Now it follows from (|3.118l) - (13.1201) and (I3.127P - (|3.128l) that Vn /2 satisfies the quadratic 
inequality: 

V n - C(l + \t\fn- l ' 2 V^ 2 - C(l + \t\fn~ 1 < 0, 



implying (|3.9p and then f|3 . 5[) . 

To finish the proof of (i) it remains to show that 

lim E{Z7 7 - 7 (i)} = v(t). (3.129) 

MOO 

In the GOE case we have by the orthogonal invariance of GOE probability measure and ()3. 105|) : 



E{Ujj(t)} = Bin' 1 Tr U{t)} -> v(t) as n — > oo. 
Besides, it follows from (13.1131) . (|3.115|) . (f2~T6l) . and (|2~20|) with p = 3, that 

„1 r 3 



(3.130) 



E{^(i)-^(t)} 



[E^ n) 



£3 



where T p are given by (|3,119p with <3?/ m = (Uji * U m j)(t, s) (cf f)3. 116j) ) . and 



5/6 n 

N < E SU P I Aln*fc»| < C(l + |t|) 5 n-^ 



A similar but much simpler argument leading to (13.127P - (13.128P allows to conclude that \T^\ < 
0(n _1//2 ), n — > oo, p = 2, 3. Hence, 

E{U j:j (t) - Ujj(t)} = 0(n- 1/2 ), n -> oo. 

This and (I3.130p yield (13.1291) and finish the proof of (i). 

(ii) The proof of (I3.107P repeats with natural modifications the one of the first part of (i). 
Namely, similarly to f|3.112|) - (|3.118|) we have for V n = Var{v n (t)}: 



V n = B{[v n (t)-v n (t)]v° n (-t)} 



it 
2 

it 
2 



f \t= 3 E nw% ] ^} - 1 3 E 



i , 4 



+ £4 



p=2 



(3.131) 



where now 



rp(n) 



pin 



1 n 

lr7 (3T^ E S+i,^ E {^^} ; P = 2 > 3 > 4 , 



,fc=l 



= t/*(t, s)v° n (-t), \D l ik (s)$ ik \ < C(\t\) 
C(\t\) is a polynomial in \t\ with positive coefficients, and 

C 4 w 6 



N < 



E 



sup 



<C{\t\)n- 2 . 



(3.132) 
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Using the argument leading to (|3.127p and (|3. 128|) it can be shown that 

2f> < C(\t\){n- l V^ 2 + n- 2 ), p = 2,3, T 4 (n) < C(\t\^-- 
This and (|3.13ip - ([3.132p allow us to write the inequality 

V n <C{\t\){n- 1 Vy 2 + n~% 



valid for any real t and implying (|3.107p . 

(iii) Statement (iii) was proved in Lemma 3.1 of |7J. 

(iv) Let t = (ti,t 2 ) and 

n 

v n i(t) = n- 1 ' 2 Ujk{ti)U kk (t 2 ) 
with U of (|3.111l) . We have similar to (|3.115|) 

v n i(t) -% 



k=l 



f 1 d n 

ll (t)= / n- 1/2 —y2u jk (t 1 ,s)U kk {t 2 ,s)ds 

■ r-l n 

= \ \ n- 1 ' 2 £ (s~ 1/2 Wff - (1 - s)- 1 / 2 W pq )^ pq (t,s)ds, (3.133) 
z Jo 1 



where now 



s) = n 1/2 5^(C^ P * U qk )(tx,s)U kk (t 2 , s) + U jk (tx,s)(U kp * U qk )(t 2 , s), 

k=l 

\D pq <S> pq (t,s)\<C(t), i€N. 
Applying ([2TT6]) and ([2720]) with p = 3 and $ = we get 



EKi(t)-w„i(t)} 



£ E{D 2 pq $ pq Ct,s)} 



2n 3/2 

p,<j=l 
n 

p,g=l 



(3.134) 



-V^ + OfrT 1 / 2 ) 



n —> oo, 



where K4 is defined in (|3.12p . Using the argument leading to (|3. 127j) and (|3. 128|) and based on 
(I2.10P - (I2.12p . it can be shown that both the terms in the square brackets of (I3.134P are of the 
order 0(n -1 / 2 ), n — >■ 00. Hence, 

EKi(t)}-E{u„i(t)} = 0(n- 1 / 2 ), ra^oo. (3.135) 

Moreover, replacing $> pq (t,s) with ^ pq (t, s)v^ ll (t) we can also obtain 

Var{v nl (t)} = E{(v n i(t) - v nl (t))v° nl (t)} = O^ 1 / 2 ), n -> 00. (3.136) 

Now it follows from (|3.135p - (|3.136p that to finish the proof of (iii) it remains to show that 

lim B{v nl (t)} = 0. (3.137) 

n— >oo 

Indeed, since by the orthogonal invariance of the GOE probability measure we have 

E{U jk (t)} = 6 jk -E{v n (t)}, 
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where v n {t) = Tr U(t), \v n (t)\ < 1, then in view of fl32) and (j2~TT|) 

ERi(ti,t 2 )} = n _1/2 E{u n (ti)}E{t; n (i2)} (3.138) 

n 

+ n" 1 / 2 ^E{c/ ifc (t 1 )C/ fc fc (t 2 )} =0(n- 1 / 2 ), n ^ oo, 



fe=l 

and we get (iv). 

(v) The scheme of the proof of (iv) is the same as the one of (iii). Namely, we have similar to 
(PHMD - Q5ZEMD (see also (ETTI5|| ): 

where now 

^ Pgm (t, «) = y^(^p * u qk )(t m , s) yi u jk (t m ',s), 

k=l m'^m 

\D l pq % qm (t,s)\<C(t), leN. 

Applying (I2.16P and (12.20P with p = 3 one can get an analog of fj3. 134H and then show that (cf 
(13T35]) - (151361) ) 

EK 2 (t)} - E{? n2 (t)} = O^ 1 / 2 ), rw oo, 

and 

Var{v n2 (t)} =0(n~ 1/2 ), n oo. 

So, it remains to prove that 

i 

lim B{v nl (t)} = TT u(i m ). (3.139) 

n— >oc A ■*■ 

m=l 

Applying (f2TTU|) and (^T2|) and then (f3T23"J) and (|2TTT]I , we get 

/ 



— E{v n2 (t)} = -±= v{w jP Uk P (ti) J] U jk {tm)) (3.140) 

fc,p=l m=2 

2 n £ / 

= E E {(*i^'fc(*i) + ( n ^™ * UjkXh)) H U jk (t m ) + E II M*™') 



n 

fc=l m=2 rra=2 m'^m 



I (Ujj(t m - si)U k k(h + si) + U jk (t m - si)Ujk(ti + si)^dsi|, 



so that by (|3,23p and ([2.1ip <9E{ti n 2(i)}/<9ii = ^(1), n — > oo, and by the symmetry 

d 



dt r 



-E{v n2 (t)} = O(l), n-^oo, m = l,..J. 



Hence, there exists a subsequence {E{v ni2 (t)}} that converges uniformly on any compact set of M. 1 . 
Now, applying the Duhamel formula (|2.8p and then (|3.140p . we obtain 



E{v n2 (t)} =E{n Ujjitm)} + -±= E E{%f/ fcp (s) J] ^ ifc (i m )}^ 

m=2 J ° ^ n k,p=l m=2 

= e| Y[ Ujj{t m )} - w 2 j ds J E{v n (s - s 1 )}E{v n2 (si,t 2 ,...,t l )}ds 1 +r n 



m=2 
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where 

,,2 rtl n , 1 



r n=~— Yl E { ( S M S ) + ( n < * U jk )(s)) H U jk (tm) 
U J ° k=l m=2 

)Ukk(s + Si) + Ujk{t m - S\)Ujk{s + S\))dsi \ 



m=2 m'j^m 

and by the Schwarz inequality, (I2.1ip . and (13. 107)) 

r n = 0(n _1 ), n oo. 
Taking into account (I3.106P we obtain that every limit of converging subsequence 

v 2 (t) = lim v ni2 (t) 

satisfies the equation: 

'o 



v 2 (t) + w 2 ds v(s - sx)v 2 (s,t 2 , ...,ti)dsx = TT v(t m ). (3.141) 

m=2 



Applying the generalized Fourier transform with respect to the variable t\ (see Proposition 12, ip , we 

get 



I 

v 2 (z, t 2 , ...,ti)(l + w 2 z~ 1 v(z)) = z' 1 Y[ v(t mJ 

m=2 

with v of (|3.4ip . Hence, v 2 (z,t 2 , ... ,ti) = v(z) \[ l m=2 v(t m ), and 

l 

v 2 (t) = Y[ v(t m ). 

m=l 

This completes the proof of (v) and the proof of the Lemma. ■ 

After this paper was completed we became aware of paper [llj by Soshnikov et al in which 
Theorem 13.41 was proved by another method and under weaker conditions. 
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